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Abstract 

The universal validity of the second law of thermodynamics is widely attributed 
to a finely tuned initial condition of the universe. This creates a problem: why 
is the universe atypical? We suggest that the problem is an artefact created by 
inappropriate transfer of the traditional concept of entropy to the whole universe. 
Use of what we call the relational IV-body problem as a model indicates the need to 
employ two distinct entropy-type concepts to describe the universe. One, which we 
call entaxy, is novel. It is scale-invariant and decreases as the observable universe 
evolves. The other is the algebraic sum of the dimensionful entropies of branch 
systems (isolated subsystems of the universe). This conventional additive entropy 
increases. In our model, the decrease of entaxy is fundamental and makes possible 
the emergence of branch systems and their increasing entropy. We have previously 
shown that all solutions of our model divide into two halves at a unique ‘Janus 
point’ of maximum disorder. This constitutes a common past for two futures 
each with its own gravitational arrow of time. We now show that these arrows 
are expressed through the formation of branch systems within which conventional 
entropy increases. On either side of the Janus point, this increase is in the same 
direction in every branch system. We also show that it is only possible to specify 
unbiased solution-determining data at the Janus point. Special properties of these 
‘mid-point data’ make it possible to develop a rational theory of the typicality of 
universes whose governing law, as in our model, dictates the presence of a Janus 
point in every solution. If our self-gravitating universe is governed by such a law, 
then the second law of thermodynamics is a necessary direct consequence of it and 
does not need any special initial condition. 
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1 Introduction 


Entropic arguments are very effective when applied to confined systems 1 in equilibrium. 
However, the universe is unconfined and far from equilibrium. Is it then correct to apply 
conventional entropy concepts to it? We think not. New issues come into play when 
the universe, which by definition is a unique system, is considered. 

We first argue that the universe’s scale and shape degrees of freedom must be treated 
differently. Scale degrees of freedom and time intervals are useful to characterize sub¬ 
systems, for which the rest of the universe provides a reference frame with respect to 
which length and time units can be defined. They lose meaning, however, when applied 
to the whole universe, which by definition admits nothing external to it. Any prop¬ 
erty of the whole universe, including functions used to quantify its typicality, must be 
expressed in terms of scale-invariant ratios of quantities within the universe together 
with the changes of these ratios. If this is done, the principles of spatial and temporal 
relationalism [1] are realized. 

Next, the specification of initial conditions is suspect: in many unconfined systems 
governed by time-symmetric laws all solutions divide into two qualitatively similar halves 
at a unique ‘central’ point. We call it the Janus point j (after the Roman god that 
looks in two opposite directions at once) and say that the corresponding solutions are 
Janus-point solutions. To our knowledge, the existence and potential significance of 
such solutions have never been noted in the voluminous literature on the second law of 
thermodynamics and its application to the universe. 

As an example of the difference that they can make, we showed in [2] that in self- 
gravitating Newtonian systems subject to natural restrictions all solutions have a Janus 
point 2 and contain arrows of time, expressed through the growth of a scale-invariant 
quantity which we call complexity , 3 that arise of dynamical necessity. Contrary to a 
widespread belief, this result shows that all the solutions of time-reversal symmetric 
dynamical systems can exhibit time-asymmetric behaviour for internal observers, who 
must be on one or the other side of the Janus point. The arrows of time are therefore 
not due to any special initial condition (or ‘past hypothesis’ [3]) but are rather a direct 
consequence of the law of the universe. This result suggests that previous discussions 
of the second law of thermodynamics have failed to note the one fact - dynamically 
enforced Janus-point solutions - that could be the true explanation for the universal 
validity of the second law in our universe. We here present more support for this view. 

First, we show that when a Janus point does exist it is the only point at which 
solution-determining data can be specified without bias. We call such data ‘mid-point 
data’. Moreover, the model that we used in [2] and call the relational N-body problem 
has properties which make it a good ‘toy’ universe and possesses mid-point data with 
remarkable properties. They are scale-invariant and non-redundant. In addition, they 
stand in a one-to-one relation to the points of a geometrical structure known as a 
projectivized cotangent bundle, which we denote by PT *S (the S stands for shape space , 
a key notion in our theory). We detail below the properties of PT *S and show that it 
carries a natural bounded measure. 

This fact plays a central role in our paper because it enables us to develop a math¬ 
ematically rigorous theory of the typicality of universes. The point is that, at least in 

1 We distinguish dynamically closed systems, on which no external forces act, from confined systems 
such as particles in a box. The former can expand freely, as is the case with the universe. 

2 In [2], this coinage had not yet occurred to us. 

3 It measures the clustering of the particles in the system 
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the relational iV-body problem, we know the law that governs the universe but not the 
probability with which any particular solution will be realized. Under these conditions, 
as argued in [4], the only rational choice one can make is to assume all solutions have 
equal probability of being realized. This is Laplace’s well-known principle of indiffer¬ 
ence. If this is to be applied consistently in a situation in which there is a continuum of 
possible solutions, the solution space must be endowed with a bounded measure. This 
is exactly what we possess by virtue of the natural measure on PT* S, whose points 
map one-to-one onto the solution space of the theory. As we will show, the existence of 
the Janus point, the scale invariance of the mid-point data and the bounded measure 
associated with them gives our theory three distinct advantages compared with the sit¬ 
uation in [4], where the corresponding measure is infinite and there is no Janus point. 
The three advantages lead us to an entirely novel overall picture of the universe. 

Indeed, the full development of our theory shows that, at least in our model, the 
universe must be described by two distinct entropy-type quantities. Both are based on 
Boltzmann’s great insight: entropy is a count of microstates compatible with a given 
value of a state function defined on the phase space of the theory. For the universe 
as a whole, we introduce the notion of entaxy (from the Greek for ‘towards order’). 
It is a scale-invariant count of microstates compatible with given values of the scale- 
invariant complexity. In any solution, it is maximal in the neighbourhood of the Janus 
point and decreases with increasing distance in either direction from the Janus point. 
The very fact that the entaxy is decreasing reflects the breakup of an initially (at the 
Janus point) very uniform universe into gravitationally bound clusters (branch systems) 
for which an increasing conventional (dimensionful) Boltzmann entropy is defined as a 
(microcanonical) count of microstates compatible with the (dimensionful) values of the 
energy and angular momentum of the clusters. Thus, there is an increasing entropy 
of the universe. It is obtained as the algebraic sum of the bound-cluster entropies, the 
overwhelming majority of which increase in either direction away from the Janus point. 
In the framework of our theory, the entaxy is clearly the more fundamental quantity: 
it must decrease, and thereby creates the conditions under which branch systems can 
form and increase their entropy. 

As a proof of principle, we consider as an example the 3-body problem and show 
that our theory is predictive: solutions that have more probable mid-point data have 
asymptotic behaviour which is characteristically different from that of the solutions with 
less probable mid-point data. Thus, observers within such a universe could in principle 
establish whether they belong to a typical or an atypical universe. This is a central issue 
in modern cosmology. Our approach suggests it may be amenable to resolution. 

We are also able to show that in all systems with Janus points j the scale-invariant 
asymptotic behaviour that develops with increasing distance in either direction from j 
is subject to dynamical attractors of the shape degrees of freedom. This leads to strong 
correlations reminiscent of the behaviour associated with retarded potentials in wave 
processes. This suggests that they too are a dynamical necessity and not the result of 
a special condition ‘in the past’. 

In this paper, we have advisedly restricted attention to the relational Wbody prob¬ 
lem. Its simplicity enables us to deduce all our conclusions with rigour and present the 
broad principles of our theory in a transparent manner. Of course, the real question 
is whether the theory is appropriate for the universe. As we note in our conclusions, 
this will very likely be so if the big-bang is not the beginning of time but instead a 
Janus-point that hides from our view a qualitatively similar twin universe on its other 
side in which the internally observed direction of time is the opposite to ours. 
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2 Shape, Size and Entropy 

Point particles moving inertially model the simplest possible universe. If the particles 
were small elastic spheres, they would model an ideal gas. If confined and in equilibrium , 
its entropy would be well defined as would changes of it under controlled adiabatic 
manipulation. But if unconfined, what can be said about the system in itself? If r™ is 
particle a’s position relative to the centre of mass R cm , the moment of inertia is 

N i N 

4 m = ' r a m = ^m a m b r^ b , M=^m a , r ab = ||r 6 - r a ||, (1) 

a a<b a 

and the root-mean-square length 


Cns ^ 1 /?. y/lgn — ^ m a rtlh r ab (2) 

y a<b 

is the degree of freedom (dof) that measures the system’s size. An external frame defines 
position and orientation dofs. Truly intrinsic to the system are its shape dofs. 

Consider now a system of inertial particles: N uninteracting particles, corresponding 
to an ideal gas in empty space. Whatever initial positions and momenta we attribute to 
the N particles, they will fly off undisturbed each along its own trajectory, dispersing to 
infinity. The moment of inertia and £ rms will grow unboundedly large. The asymptotic 
configuration will have £ ims = oo and a shape determined only by the momenta. The 
information about the initial positions will be completely erased, and the cm coordinates 
r™ and momenta p“ m will become more and more correlated as £ Iins — >■ oo: each particle 
asymptotes to recession from R cm ; the angles Z {r™, p“ m } — > 0 and any two particles 
will tend to separate with speed v ab oc r ab as in Hubble-type cosmological expansion. 

We now note that our system of inertial particles has time-reversal symmetric equa¬ 
tions of motion: if we evolve those same initial conditions backwards in time (or, equiv¬ 
alently, reverse all the initial velocities), we obtain the same asymptotics at t —> —oo, 




-t 

Figure 1: A planar solution of an 8-body system of inertial particles. Each particle moves 
at a constant speed in a straight line (from the perspective of an inertial reference frame and 
with Newtonian time [1]), r a = r a (to) + v a (i — to). At both time boundaries the system tends 
to two parity-conjugate shapes, which are completely determined by the velocity vectors v a 
and do not depend on the initial (or better mid-point, see below) coordinates. Because the 
system is in two dimensions and not three, the parity conjugation actually leaves the shapes 
identical. 
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with the difference that the asymptotic shape is the parity-reversed image of the other. 
Moreover, £ rms tends to oo both at t —> oo and t —>■ —oo. It is easy to see that f r 2 ms is a 
positive quadratic form in t, and therefore reaches a unique minimum in each solution. 
For this reason, systems of inertial particles are the simplest with a Janus point j: £ ims 
always passes through a unique minimum, and £ rum —» oo as t —> ±oo. It would clearly 
be inappropriate to try to specify ‘initial’ conditions freely at either t —> ±oo limit due 
to the tight correlations between r™ and p“ m . 

Entropy growth is often described qualitatively as passage from an atypical (ordered) 
to a typical (disorded) state. Let us consider this in the light of the observations we 
have just made. If we specify uncorrelated initial data for a system of inertial particles, 
the state corresponding to those data must be near the Janus point of the solution they 
generate - at any distance from j the particle motions become correlated, as we have 
just noted. Thus, the state of the system of inertial particles becomes less typical going 
away from j, and its ‘arrows of typicality growth’ point towards j. This is typical of all 
unconfined systems with Janus-point solutions. 

It is illuminating to compare this behaviour with that of an ideal gas in a confined 
space. In contrast to the system of unconfined inertial particles, the accessible phase 
space of the system then has a bounded measure and is subject to Poincare recurrences: 
its phase-space point will return arbitrarily often arbitrarily close to any point previously 
visited. Although the time between close returns is extremely long (for a mole of gas 
much longer than the age of the universe), this behaviour is inevitable due to the fact 
that the accessible phase space is bounded. Because of the recurrences, the (Boltzmann) 
entropy S B of the system must fluctuate. For most of the time, S B will remain very close 
to its maximum, which corresponds to thermal equilibrium (heat death). However, there 
will inevitably be occasional fluctuations to lower values of S B . These will be ever rarer, 
the deeper the dips in S B . 

The plot of 5 b against the time t will exhibit qualitative time-reversal symmetry, 
reflecting the time-reversal symmetry of the law that governs the system. In a famous 
letter to Nature [5], Boltzmann conjectured that we find ourselves in a world in which 
entropy increases because we happen to be on one side of a deep Poincare dip and define 
the direction of time as the direction away from the dip, towards equilibrium and heat 
death. Intelligent beings on the other side of the dip will find that their time flows in 
the opposite direction, again to equilibrium and heat death. Since typicality increases 
with entropy, in such a system the arrows of typicality growth point away from the dip. 
Because the behaviour is symmetric about the dip, there is qualitative time-reversal 
symmetry about it just as there is for the unconfined system of inertial particles. The 
crucial difference is, however, that the arrows of typicality growth in the confined system 
point away from the dip, and not toward the dip as the Janus-point system. Moreover, 
whereas there is a single Janus point in the unconfined system, there are infinitely many 
dips in the confined system. 

The origin of the ‘towards-away’ mismatch and also the fact that in one case there 
are inhnitely-many dips but in the other just one special point is obvious. If a system of 
elastic inertial particles with £ rms increasing freely is abruptly confined in an elastic box, 
collisions of the particles with the box and each other will rapidly reverse its atypicality 
growth, and a statistically typical equilibrium state with Poincare recurrences will be 
established. Prevention of £ ims growth forces the shape dofs to equilibrate and explains 
both of the mismatches. 

Natural processes whose time reverse is never observed are said to be irreversible 
and associated with entropy growth. An example sometimes given ([6], p. 33) is that 
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“a gas will explode into a vacuum, but will never spontaneously implode into a smaller 
volume”. This is certainly true within the universe, but a universe made only of N 
inertial particles ‘implodes’ from t = —oo to j and then ‘explodes’ to t — +oo for either 
nominal time direction. 1 

Now consider the role size plays in ideal-gas entropy. Up to an additive constant 

S = Nk log (UT 3 ^ 2 ), (3) 

where k is Boltzmann’s constant, V the volume, and T the temperature. Under re¬ 
versible adiabatic expansion, 5 S increases because the system constantly re-equilibrates, 
thereby maximizing the entropy subject to new constraints. Note that, the particles and 
container being elastic, T remains constant; S increases because V does. But without 
the constraint-enforced shape equilibration, calculating an entropy is problematic. The 
Gibbs H function is constant under Hamiltonian evolution: it is extremalization of H 
subject to constraints that leads to experimentally verifiable entropy values [7]. The 
Boltzmann entropy S B , based on a one-particle distribution in a six-dimensional phase 
space and obtained from the Gibbs H by integrating over the degrees of freedom of all 
but one of the particles, can increase. However, it inherits the ensemble nature of the 
Gibbs H function, and it is ensemble spreading that leads to the appearance of typicality 
growth. 6 In reality, every single microsolution within the ensemble tends with increas¬ 
ing distance from j to an ever more atypical, increasingly correlated state. This is quite 
unlike a confined gas as it approaches equilibrium, in which virtually all correlations are 
being destroyed in every microsolution. 

Moreover, phase-space volume is a product of size and shape factors. In unconfined 
Janus-type systems, the sizes in any Gibbs ensemble can spread freely, so the size factor 
will grow in both time directions. To keep H constant, the shape factor must decrease. 
The shape dofs are therefore subject to dynamical attractors , as demonstrated clearly 
in Fig. 2. Their existence is hidden unless scale is factored out. 

Finally, overall size, unlike dimensionless shape dofs, only has meaning relative to an 
external scale, which does not exist for a universe. The justification for saying that our 
universe is expanding is based on ratios, in the first place of wavelengths in red shifts. 
Observations of them make it meaningful to say that the universe at a given epoch is 
twice as large as it was at some earlier epoch, but not to ascribe size to any specific 
epoch. Relative scale is meaningful but absolute scale is not. 

In the light of these considerations, of which we have found little recognition in 
the literature, 7 we conclude that entropic arguments require revision when applied to 
the universe. We can summarize our position as follows: in any dynamically closed 
but unconfined system that can increase arbitrarily in size (in the sense defined above) 
and whose solutions have a Janus point, judgements as to the growth of typicality 

4 A half-solution, from Urns = oo to j, has a seemingly implausible initial condition at £ rms = oo 
like an advanced-wave solution in electrodynamics. This suggests that the invariable observation of 
retarded waves in the universe is due to its having a Janus point j and our observations being made in 
an asymptotic region. 

5 Adiabatic expansion of a gas in a piston is reversible because the overall entropy of the sys- 
tem+environment does not change, even though the gas increases its entropy in the process. 

6 Note also that the Boltzmann entropy does not in general (unless one is considering an ideal gas) 
lead to the correct empirically determined equilibrium thermodynamic relations [7] because its one- 
particle distribution function washes out correlations taken into account in the Gibbs H function based 
on the N-particle distribution function. 

7 Sloan [8] points out that in cosmology the presence of the scale factor in the case of a minimally 
coupled scalar field hides the presence a shape attractor. 
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Figure 2: This diagram shows a bundle of 30 solutions of the inertial planar 3-body problem 
in the space of scale-invariant configurational degrees of freedom, the shape sphere (defined in 
the next section and explicitly constructed in Appendix A). Planar means that the angular 
momentum is perpendicular to the plane of the 3 bodies, and therefore the motion unfolds 
on a fixed plane. Each solution is divided in two (a red and a yellow part), at each side 
of its Janus point. All solutions have the same momenta but the mid-point coordinates are 
normally distributed around a central shape, with a variance of ~25%. In this way one gets 
a distribution in shape space whose evolution, at both asymptotic ends, converge to a single 
shape (as noted in the caption to Fig. 1, in two dimensions parity-conjugate triangles are 
identical), completely determined by the momenta. The presence of dynamical attractors is 
here manifest; later, when we add gravitational interaction to inertial motion, we will find an 
even more impressive manifestation of such attractors. 

need to be based on the behaviour of the shape dofs, which are the only ones that are 
directly observable. The ‘invisible’ (and single) scale dof does play an important role in 
‘allowing’ the shape dofs to change their overall typicality but should not be included 
in any calculation of typicality. 

Thus, the critical issue in the definition of entropy for the universe is, as we see it, 
this: as a physical system, is the universe confined or unconfined? 8 If the latter, which 
seems intuitvely natural, how can we define for it a bounded measure, without which 
an entropy-type count of microstates for a given macrostate is ill defined? 

3 The Relational iV-Body Universe 

We obtain a model of the universe that is much more realistic than inertial particles 
by introducing Newtonian gravitational interaction between particles of mass m a , a = 
1... N. It is a model ideally suited for our purposes: simple but nontrivial enough to 
illustrate our main points. It is also valuable in that it brings out the difference between 
gravitational and nongravitational forces. 

8 Let us stress that our notion of confined is not equivalent to finite, or compact. These are two 
completely unrelated concepts. 
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What we mean by a relational model is one that makes no use of extrinsic elements. 
The dimensionless mass ratios m a /M and the dimensionless inter-particle separations 
s a b '■= f' ab / L, are intrinsic quantities. If N > 5 these are not all independent; the 
instantaneous shape is determined by 3 N — 7 independent combinations of the s a b. The 
extrinsic elements in Newtonian dynamics enter through an inertial frame of reference. 
This provides an external time parameter t, a centre-of-mass velocity, 9 an orientation of 
axes and a scale. None of the extrinsic elements are directly observable. The 3 N particle 
coordinates r a in an inertial frame are therefore heterogeneous, mixing instrinsic and 
extrinsic degrees of freedom: 

• three locate the position of the system in the frame, e.g., the centre of mass 

N N 

Rcm = m a r a/M, M = ^2,m a (4) 

a= 1 a= 1 


• three fix its orientation, e.g., the orientation 10 of the 3 eigenvectors of the centre- 
of-mass moment of inertia tensor (fixing the ‘principal axes of inertia’) * 11 

N 

l = XI r r ■ r r - r r ® o, c = ^ - r cm , (5) 

a= 1 


one hxes its size, measured for example by the moment of inertia (1): 

1 N 1 

I cm == ^m a r c a m -r c a m = —^2m a m b r 2 ab , r ab = ||r 6 - r 0 ||, (6) 

a =1 a<b 


• the remaining 3 N — 7 are the relational, scale invariant data [9, 10, 1] (in short, 
shape data). They are angles and ratios of the separations r ab between the parti¬ 
cles. Only these are observables. 

None of the extrinsic elements are directly observable. To establish how they affect 
what is observable, we define shape space S = Q/Sim by quotienting the Newtonian 
configuration space Q ~ R 3jV by the similarity group Sim: 

Sim := |r a ~ r' a iff 3 6 G R 3 , 0 G 50(3), ip > 0 / r' a = <^(f2 r a + 6) } . (7) 

Any Newtonian solution can be projected to a curve in S. This is simply a succession 
of shapes without the earlier-later ordering provided by t. In the 3-body problem, the 
simplest non-trivial W-body problem, two angles, a and /3, determine the instantaneous 
shape of the triangle that the particles form. The direction of change is da/d/3. In a 
purely intrinsic dynamics, a, f3, da/d/3 should determine the evolution a = a(/3). This is 
not so for a generic Newtonian solution because five dimensionless quantities of extrinsic 
origin play a role. They are: 

r> By Galilean relativity, this cannot be observed in the behaviour of the s ab and is of no concern. 

10 Notice that the orientation of a vector involves just two degrees of freedom, but the orientation of 
a system of three orthogonal vectors involves three dofs. 

11 This is just an illustrative example. It is well known that there are spherically- or cylindrically- 
symmetric configurations whose principal axes of inertia are degenerate and cannot be used to fix an 
orientation in space. However, these configurations are special (measure zero). 
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• the direction of the cm angular momentum vector L = ^ a r[j m x p“ n (two dofs), 

• the ratio T cm /V Nevi of the cm kinetic energy to the Newtonian potential energy, 

• the ratio T R /T S between the rotational and shape contributions to T cm , 12 

• the ratio T D /T S between the dilatational and shape contributions to T cm . 

The values of the five dimensionless quantities are not encoded in the shape s and 
direction d of the curve in S but in the way the curve bends at s. Intrinsic s, d initial 
data fail to determine the intrinsic evolution (they fall short by five for all N > 3). 

The three data related to the angular momentum L (two angles and the ratio T R /T S ) 
are present because orientation is absolute in Newtonian theory. They cease to play a 
role if L = 0. Similarly, if the total energy E vanishes, E = 0, the absolute time plays 
no role and the ratio T cm /V New must have the value —1. Otherwise, it will change along 
the curve in S and can have any value at an initial shape s. It is important that the 
conditions E = L = 0 are conserved by Newtonian dynamics. This means we can 
consistently eliminate the effect of four of the five extrinsic quantities. 

We now come to the most critical point: whereas the angular momentum L is con¬ 
served in the iV-body problem, the dilational momentum 13 D = ^E a r a m ' Pcm is not. 
(Unlike L, D does not commute with the Hamiltonian.) Because of this, the ratio 
T d /T s of the dilatational and shape kinetic energies cannot vanish throughout the evo¬ 
lution; even when E — L = 0, a shape and direction in S are in general insufficient to 
determine the evolution curve in S. What one finds is that for any point s in S and 
any direction d at s there is a one-parameter family of projected Newtonian solutions 
through s whose tangents all have the same direction d. With increasing distance from 
s, they ‘peel off’ from the common direction d into manifestly distinct solutions. This 
characteristic behaviour is due to the role that relative (as opposed to absolute) scale 
plays in Newtonian gravity, as we shall now show. 

For any potential V homogeneous of degree k, 14 it is easy to show that 

I cm = 2E~(k + 2)V. (8) 

The Newton potential is homogeneous of degree —1, so that for it J cm = 2 E — F New . 
Moreover, Wew is negative-definite, so that if E > 0 

J cm > 0. (9) 

This means that 7 cm is concave upwards as a function of t and (half) its first deriva¬ 
tive, 

1 N 

2 J - = E rT-Pl = A ( 10 ) 

a =1 

is a monotonically increasing function (under time reversal D changes sign and therefore 
is increasing in both nominal time directions). 15 This is a consequence of the Lagrange- 
Jacobi relation. At one unique point j we have D = 0, the Janus point of the system. The 

12 The velocity decomposition theorem [11] states that at any instant T cm = T R + T D + T S . 

13 This name for D, which is the generator of dilatations, was introduced in [9]. It has the same 
dimensions as angular momentum, the generator of rotations. There is no name for D in the IV-body 
literature; it is usually denoted by J, probably for Jacobi. 

14 A function f(xi) is homogeneous of degree k if f{axf) = a k f(xi) for any constant a > 0. 

15 It is an immediate consequence of the existence of this monotonically inceasing (Lyapunov) function 
that the solutions of the relational N-body problem cannot be periodic. The system is not ergodic and 
cannot be subject to Poincare recurrences. 
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value of 7 cm at D = 0 is purely nominal since it depends on the unit chosen to measure 
distance. However, the ratio of the J cm values at distinct points in the evolution is 
physical and can in principle be determined from the curve in S. This is also true for 
the instantaneous value of the ratio T D /T S and the point where both it and D vanish. 
Overall (absolute) scale has no objective meaning but relative scale does. 16 

Thus, the basic structure of Newtonian evolution is not scale invariant. It is however 
remarkable that each individual, observationally distinct evolution curve in S corre¬ 
sponding to Newtonian solutions with E = L = 0 is uniquely and nonredundantly 
determined by scale-invariant quantities. This is due to the existence of the single point 
at which D — 0. As we have noted, this is a Janus point, and at this point in the 
evolution there is no dilatational kinetic energy: T D = 0. Thus, at D = 0 a point and 
direction in S do determine a solution, moreover non-redundantly. 1 ' We have mid-point 
data. 18 When working in phase space, we will call the intersection of the hypersurfaces 
E = L = P = 0 19 and D = 0 the Janus surface , which we will denote by the symbol J. 

The manner in which these mid-point data are scale invariant is worth spelling out: 
for N particles, 3N — 7 numbers fix the shape s and 3N — 8 the direction d (in any space, 
a direction is fixed by one less number than a point). The conditions E = L = 0, which 
fix the class of solutions considered appropriate for a dynamically closed universe, and 
the condition D = 0, which fixes the point at which the mid-point data are specified, are 
scale invariant. Since we regard as objectively different only those Newtonian solutions 
that map to distinct curves in shape space S, we conclude that our space of solutions 
has 6N — 15 dimensions. Finally, we pointed out earlier that through every shape s G S 
and for every direction in S at s there passes a one-parameter family of solutions all 
tangent to d. Just one of these solutions has D = 0 at s. 

The non-redundance of mid-point data, which leads to the odd-dimensionality of the 
solution space, is due to dynamical similarity , which arises because the Newton potential 
is homogeneous [14, 15]. By rescaling all the initial (or better ‘mid-point’) momenta 
p“ —> kp a by a common constant factor 7, one obtains a solution that is congruent 
(geometrically similar) to the original one, namely r a —>• k~ 2 r a . If parametrized with 
Newtonian time (which can be defined relationally as the ephemeris time of [16, 1]), this 
solution runs at a different pace, which is just a rescaling t —* k~ 3 1. So the two solutions 
differ by an unobservable global rescaling and an equally unobservable time rescaling. 

16 To create a dynamics of pure shape [9, 1], in which one can simultaneously have D = E = L = 0 
and therefore a point and direction do determine the evolution in S, the Newton potential must be 
replaced by one that is homogeneous of degree —2. 

11 As emphasized in a valuable comment to one of us (JB) by Halliwell [J Halliwell, private communi¬ 
cation, September 2009], a general phase-space point determines a solution and, redundantly, an initial 
point on it. The existence of the unique point at which D = 0 enables us to specify solution-determining 
data at that point nonredundantly. 

18 In [12] Chen and Carroll contemplated time-symmetric solutions of time-symmetric equations with 
entropy growing to infinity either side of a unique minimum. Janus-point systems confirm this ‘middle 
hypothesis’ [13] though, as we will see below, our entropy-type quantity has a maximum, not minimum, 
at j. As Carroll [13] notes, it is intuitively desirable that systems with time-reversal symmetric equations 
should have solutions that are qualitatively time symmetric. In fact, this has long been recognized to 
be so in confined systems in which infinitely many Poincare recurrences occur. It is strange that, 
in the arrow-of-time literature, the common occurrence of systems with Janus point j and just two 
qualitatively symmetric branches either side of j does not seem to have been noticed. Besides inertial 
systems and the relational TV-body problem (the TV-body problem with E = L = 0), two other obvious 
examples are the TV-body problem with replusive gravity and (liquid or electromagnetic) waves with 
compact support. 

19 P = p a = 0 is the total linear momentum of the system. 
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They are indistinguishable to observers living in such solutions, who have access only 
to ratios between space and time intervals within their universes. 

Although S aids intuition, calculations are more convenient in the Newtonian phase 
space T. This contains the following seven constraint hypersurfaces: 

N N N 

p = p“ = °> L = Ta x p a = o, D = r « ■ p “ = °’ ( n ) 

a =1 a= 1 a= 1 

According to the classification of constraints introduced by Dirac [17, 18], these seven 
constraints form a first-class system; the intersection J C T of all the mentioned hy¬ 
persurfaces is a 3 N — 7-dimensional manifold, over which the constraints define seven 
independent Hamiltonian vector fields which can be integrated into seven-dimensional 
gauge orbits. The seven constraints close the Lie algebra of the similarity group Sim = 
R + ix ISO( 3). Except at some singular points which can be ignored because they form 
a set of measure zero (according to the measure we are going to introduce below), one 
and only one gauge orbit passes through each point on J, and therefore the quotient 
space of £ by the orbits (‘reduced phase space’) T red = J/Sim is well defined [18] and is 
a symplectic manifold. In our case we can call this space shape phase space. 

The Hamiltonian constraint 





+ v„ 


E 


cm 5 


14 


£ 

a<b 


m a m b 

r a b 


( 12 ) 


is different: it generates the dynamics of the system and is first-class with respect to L 
and P. But H and D are second-class, 


{77, D} = 2 Y, + ^ew = 2 77 - P New + 0, 

ZTfl a 


(13) 


which is a quantity that never vanishes. This is the phase-space formulation of the 
Lagrange-Jacobi relation presented above. Therefore each dynamical orbit generated 
by 77 (the integral curves of the Hamiltonian vector field associated with 77) intersects 
D = 0 and J only once, at the Janus point. In the next section we will use this fact to 
specify mid-point data at the surface J, and to identify the space of solutions (and the 
measure on it) with T red . This will not be the end of the story, because the solution space 
thus defined will still be non-compact. The final step will be to use dynamical similarity 
(realized, in phase space, by the equivalence relation p a ~ k p Q ) to further quotient T red 
and obtain the space of physically distinct solutions PT *S (the notation will be clear in 
the following), which is compact with respect to a naturally defined measure on it. 


4 Defining Typicality of Microstates in Cosmology 

A review of the conditions that make it possible to define entropy and typicality in stan¬ 
dard situations in thermodynamics and statistical mechanics makes it clear that some 
modification of concepts must be made for the universe. We have seen little discussion 
of this issue in the literature except for an approach initiated by Gibbons, Hawking and 
Stewart [4] in an attempt to quantify the typicality of cosmological solutions that lead 
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to a given number of inflationary e-foldings. 20 We follow their basic procedure but with 
important modifications that enable us to avoid Schiffrin and Wald’s criticisms [19] of 
the GHS approach. 

• GHS want to define a measure on the space of solutions, assuming a Hamiltonian 
system with time-independent Hamiltonian. 

• Each solution will lie on a constant-energy surface H(p, q) = E on which there is 
a well-defined induced Liouville measure, so they fix the energy once and for all 
and define a measure on the solutions with energy E. In our case E = 0. 

• In order to count each solution once, they assume the existence of a surface that 
intersects the constraint surface H = E transversely and once only for each solu¬ 
tion. In general systems one is able to guarantee these two conditions only locally 
(in a neighbourhood of each phase-space point). We have a natural choice for this 
surface: D — 0, which is a good global gauge-fixing of H = 0 because it intersects 
it transversely and once only per solution, as we proved above. This is a very 
significant advantage not present in [4], 

• Having obtained their intersection surface, GHS proceed to prove that the 2-form 

cu n _i induced on it by the natural symplectic form of phase space co n = d Pi A d q l is 
invariant under the Hamiltonian phase flow (the Lie derivative with respect to the 
Hamiltonian vector field associated to H). This implies that the associated volume 
form on the constraint surface ^ (the (n — l)-th exterior power of oo n -i) is 

invariant too, and consequently it is independent of the choice of intersecting 
surface. 

These results can be better expressed with the language of constrained Hamiltonian 
systems: the measure is 


Hn-i = | det{y, H}\5(x)5(H)Y[dpidq k , (14) 

i,k 

where | det{y,i/}| is the Faddeev-Popov determinant (in this case, with only a pair of 
second-order constraints; it is just the absolute value of the Poisson bracket between H 
and x). The measure p n -i is invariant under the Hamiltonian phase flow. 

Moreover if there are further constraints H a , a = 1 ,...,m, as in our case, one 
specifies an identical number of gauge-fixings x.a, and the measure 

A I>n-m = I det{Xa, H b }\Y[fi{Xa)5(H b )Y[dpAq k (15) 

a,b i,k 

does not depend on the choice of gauge-fixings. 21 

20 All of these studies are of the massively symmetry-reduced FLRW models with a homogeneous 
scalar field. All shape (conformal) degrees of freedom are frozen. In contrast, we aim to develop a 
theory of the typicality of universes without any symmetry reduction; our only restriction is that the 
universe is assumed to be spatially closed (and globally hyperbolic). 

21 The constraints and gauge-fixings = (Hi,..., H m ,xi, ■ ■ ■ ,Xm), must satisfy some regu¬ 
larity conditions (the Jacobian matrix must have rank 2m, or alternatively the gradients 

dHi, ..., dH m , dyi)..., dXm must be locally linearly independent on the constraint surface, so that 
dHi A • • • A d H m A dyi A • • • A d Xm ^ 0) (see [18] sec. 1.1.2). 
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The last remark allows us to improve on the GHS measure. In fact the second-class 
pair 

= + o = E r «p“- < 16 > 

a a 

has the special property that each of its two members taken separately closes two distinct 
first-class systems, (H, P,L) and (_D,P.L), with the spatially-relational constraints P 
and L: 

{H. P} = 0, {H, L} = 0, {L\ P>} = e ijk P k , (17) 

{L\ Pi} = e i i k P k , {P, D} = P. {D, L} = 0. (18) 

This in turn implies that the measure 

a H = \det{(D, X p,XL),(H,P.L)}\5(D)S(x P )S(XL)S(H)S(P)S(L)]Jdp i dq k , (19) 

i,k 

which is based on gauge-fixing of the system (H, P. L) by D (xp and xl are appropriate 
gauge-fixings of P and L J2 ), is identical to the measure 

^ = \det{(H, XP ,Xi,)AD,P:m^H)S(xi>)S(XL)S(D)5(P)S(L)l[dp i dq k , (20) 

i,k 

based on gauge-fixing of the system (D, P. L) by H. In fact, on the constraint surface 
P = L = 0 we can write the first Faddeev-Popov matrix as 

( {H,D} {H, X p} {H, X l}\ ( {H,D} {H,x p} {H, X L } \ 

{P ,D} {P , Xp} {P , Xl} « 0 {P • Xp} {P,Xl} (21) 

V {L, D} {L, xp} {L, Xl} / \ 0 { L -Xp} {L,Xl} J 

because of course {P.H} « {L, H} « 0, and the other Faddeev-Popov matrix as 

( {D,H} {D, X p} {D, X l}\ / {D,H} {D, X p} {D, X l} \ 

{P,^} {P,Xp} {P) Xl} = 0 {P ■ Xp} {P-Xl} (22) 

V {L ,H} {L, xp} {L,x l } / \ 0 ( L ^p} {L, xl}/ 

(in this case because {P,H} ~ {L,H} « 0), and their determinants will coincide in 
absolute value: 


/ { D,H } 

{Axp} 

{D, x l} \ 

( {H, D} 

{H, Xp} 

{H, xl} \ 


0 

{P,Xp} 

{P- Xl} 

= — det I 0 

{P-Xp} 

{P, Xl} 

. (23) 

V 0 

{L-Xp} 

{L- Xl) / 

V o 

{L.Xp} 

{L-Xl} / 



22 A very natural choice of Xp and Xl is the following: xp = R- cm = jj J2 a m a r a, the coordinates of 
the centre-of-mass, and xl = (1™^™^™)' where 

N 

JU = E (1 r r • rr - rr ® rD , 

a—1 

is the inertia tensor. This last gauge-fixing requires the inertia tensor to be diagonal, so that the 
axes of the system are aligned with the principal axes of inertia. This gauge-fixing is not always well- 
defined: in collinear configurations, for example, it fails because the inertia tensor has rank 1, but these 
configurations are of measure zero and need not concern us here. The Faddeev-Popov determinant of 
the second-class system (P,L,xp,Xl) is just (proportional to) detl cm , the determinant of the inertia 
tensor. Then xp and xl are first-class wrt both D and H (the only nontrivial one is {H, xl}, which 
is a linear combination of D and L). 
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We see that the fact that the two constraints D and H separately close first-class systems 
with the gauge constraints P and L is crucial for the equivalence of the two measures. 
This property is what lies behind the mechanism called symmetry trading , which is at 
the basis of Shape Dynamics [20, 1], and we see it here realized. It gives us a further 
significant advantage compared with [4]: direct access to a gauge-invariant measure. 

The system (H, P,L) now describes the solutions of the E = L = 0 Newtonian 
iV-body problem in a temporally relational fashion, 23 and the gauge orbits of (.H , P. L) 
on the constraint surface H = P = L = 0 are solutions of the relational IV-body 
problem. On the other hand, the system (Z4,P.L) is a purely kinematical description 
of IV-body shapes and their conjugate momenta; 24 it knows nothing of dynamics and 
solution spaces. Following Arnold [21] (Appendix 5), we notice that the system (D. P. L) 
realizes a Poisson action of the Lie group Sim on the A r -body phase space T. Then there 
is a unique construction called reduced phase space which gives a synrplectic manifold 
T red obtained by quotienting the constraint surface Z4 = P = L = 0by the gauge orbits. 
Our measure op, is just the volume form associated with the synrplectic form of T red . 
Moreover, Arnold proves (a slightly more general version of) the following theorem: 

The reduced phase space T red is symplectic and diffeomorphic to the cotangent 
bundle of the quotient configuration manifold S (shape space). 

In particular, the canonical symplectic form of the cotangent bundle is equivalent (sym¬ 
plectic) to that of reduced phase space, and therefore the volume form ap is given by 
the canonical one on T* S. This construction is completely natural and unique. 

To summarize, we have proven the equivalence of our GHS-type measure on the 
space of solutions of the E = L = 0 Wbody problem to the measure on the reduced 
phase space T red of the scale- and spatially-relational system (74. P . L ), which we call 
shape phase space. Moreover, this measure is the same as the canonical volume form on 
T* S, the cotangent bundle to shape space. 

But we still do not have exactly what we need: the measure on T* S is infinite. In fact 
one can easily see that shape space S is compact (gauge-fixing D by fixing the moment 
of inertia results in a hypersphere), but its cotangent bundle is not compact according 
to the canonical measure: one can take arbitrarily large cotangent vectors. However, as 
we pointed out above, there is a further redundancy: rescaling the shape momenta just 
gives a similar solution (meaning rescaled by the same amount at all times). It is then 
clear that, if we mod out by momentum rescalings (and reflections, because changing 
the sign of the mid-point momenta gives the same solution, just pointing towards the 
opposite time direction), we will end up with a compact space. We can call this the 
solution space (or the space of physically distinct solutions). It coincides with PT* S, 
the projectized cotangent bundle to shape space. Thus, although the relational Wbody 
problem is not fully scale-invariant, we have a complete non-redundant scale-invariant 
representation of its solution space. 

However, PT* S is not a symplectic manifold, because it is odd-dimensional. As the 
final step in establishing the mathematical structures we need, we would like to define 
a measure on it, but the measure on T* S does not project uniquely to a measure on 
PT* S. Some more input is needed. Here there must be a choice on our part, because 
the natural geometrical structures possessed by the various configuration spaces we 
considered so far do not suffice to specify uniquely a measure on PT* S. One way we 

23 The Hamiltonian being a constraint embodies a description that does not depend on any particular 
time parametrization. The system (H, P,L) realizes spatial and temporal relationalism [1]. 

24 The system (D,P,L) in its turn embodies scale and spatial relationalism [1]. 
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could induce a measure on PT *S is to fix the norm of cotangent vectors to 1. This gives 
the sphere bundle 7j*S [22], which is the double cover of PT *S (each point in PT *S maps 
to two points related by momentum reflection in T* S). But in order to have a norm on 
cotangent vectors we need a metric on S, and this is where our choice enters. 

The extended configuration space Q is endowed with a natural metric, the mass 
metric [1], 

N 

ds 2 = ^2 m a dr a • dr a , (24) 

a= 1 


and, as shown by Montgomery [23] (p. 320), this induces a metric on the quotient of 
Q by the isometry group of ds 2 (a construction called Riemannian quotient). But the 
metric above is not invariant under dilatations, and therefore the Riemannian quotient 
procedure does not apply to S = Q/Sim. We need a metric whose isometry group is Sim. 
Any metric that is conformally related to (24) through a factor that is homogeneous of 
degree —2 in r a would suffice. In this paper we appeal to the principle of simplicity and 
make the choice 


ds 2 = 


N 

£ 

a= 1 


dr„ • dr,. 


m n 


(25) 


which is the simplest we can imagine. Any other choice would seem to lack a sufficient 
reason. Ultimately, our choice, for all its simplicity, must (like any other choice) be 
tested through the predictions it makes for observable quantities. We will discuss this 
below. 

The metric (25) gives the following norm for a cotangent vector p a : 


N 


ion y ' 

a =1 


p a p a 

m a 


(26) 


which we may set equal to one, obtaining a measure on Tj*S and consequently on PT* S. 
This measure is finite and probabilities are well-defined with it. 

For illustration, in Appendix A we flesh out the 3-body case in all its details, showing 
explicitly the analytic form of the measure on PT* S. 


How we avoid Schiffrin and Wald’s criticisms 

In the remainder of this section, we indicate how we avoid the two most serious difficulties 
identified by Schriffrin and Wald [19] in the approach initiated by Gibbons, Hawking 
and Stewart [4], 

The first of these relates to the conditions under which a count of microstates is 
physically justified and well defined. If the microstates are discrete and finite in number, 
the count of them is unproblematic. However, if there is a continuum of microstates, 
the count must be based on a measure, which must be finite if unambiguous statements 
are to be made. This is an essential condition. Schiffrin and Wald point out that the 
various studies based on the approach of [4] gave discordant results because they used an 
infinite measure, in which case probabilities become hopelessly ambiguous. 25 At least 
in the A-body problem, we resolve this problem because we work on the solution space 

25 The probability of a finite-measure region would always be zero. That of a region whose complement 
has finite measure would always be 1, and an infinite-measure region whose complement has infinite 
measure too would have an indefinite probability. 
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PT* S, which has a finite number of degrees of freedom and, as we have just shown, can 
be ascribed a finite measure. 26 , 2 ' In the field-theoretic context of general relativity, in 
which the classical theory has infinitely many degrees of freedom, it may be necessary 
to assume [19] that, in some as yet unknown way, quantum gravity supplies an effective 
cutoff. 

The final conceptual issue that we must discuss concerns our basic statistical as¬ 
sumption concerning the probability with which particular solutions are likely to be 
realized. The situation is this. We assume we know the law of the universe and that, as 
in the IV-body problem with E = L = 0, all of its solutions have a Janus point. Then 
each point on PT *S stands in a one-to-one relation with one of the distinct solution 
curves in S. We know the irreducible miminum of data needed to fix a solution. What 
we do not know is the probability with which any particular solution will be realized. 
The law by itself gives no guidance on this matter. 

Under these circumstances there is no rational stance we can adopt other than 
Laplace’s principle of indifference: 28 a rational agent who knows a certain number 
of situations can be realized but has no other information can do no better than ascribe 
equal probabilities to the various outcomes. This was clearly good advice for any gam¬ 
blers whom Laplace might have advised on betting tactics. In fact, in Bayesian data 
analysis the principle of indifference is widely used to eliminate bias. 

We are, de facto, adopting the fundamental postulate of statistical mechanics, ac¬ 
cording to which the occupation of any state is assumed to be equally probable to that 
of any other state, i.e., Pi = 1 /fi, where is the number of accessible states when 
their number is finite. One is usually justified in making this assumption for an isolated 
system for which one has good reason to believe that it is in equilibrium. One assumes 
that the system is ergodic, so that averages over time can be replaced by averages over 
phase space. It is obvious that we are applying the postulate under entirely different 
circumstances. 

Indeed, one might think that nothing can be predicted on the basis of ignorance. 
However we will show that, applied appropriately, it has the potential to lead to strong 
predictions for the mid-point data, the consequences of which can in principle be tested 
observationally. We say fin principle’ because the actual calculations needed to arrive at 
definite predictions for observable quantities are subject to uncertainties, above all ones 

26 There are two aspects of our measure that require a good motivation. One, our choice of (25), we 
have already discussed. The other is our use of the Liouville measure induced on T* S. This appears 
as a natural, indeed obvious choice. The very possibility of defining Newtonian dynamics of point 
particles interacting through universal gravitation derives ultimately from the geometrical properties of 
Euclidean space and the notion of mass. Taken together in the form of the mass-weighted coordinates 
sjm a r a , they lead directly to the kinetic metric m a dr a ■ dr a that underlies Newtonian dynamics. 
The rules for variation of the action lead to the distinguished status of the canonical momenta and 
the definition of the symplectic form. Thus, the Liouville measure directly reflects the fundamental 
structure of phase space, which (as we have noted) itself derives from geometry and the notion of mass. 

2 ‘ In loop quantum cosmology a ‘bounce’ replaces the big-bang singularity of classical general rela¬ 
tivity. Ashtekar and Sloan [24] exploit the measure induced at the bounce to specify initial conditions 
at the ‘bounce’. This enables them to define a measure on the space of solutions of symmetry-reduced 
FLRW solutions with a homogeneous scalar field and estimate probabilities for the number of inflation¬ 
ary e-foldings. The bounce is clearly a quantum Janus point. In contrast to [24], we restrict ourselves 
in this paper to purely classical theories with a Janus point and do not make any symmetry reduction, 
which in the case of the N -body problem would correspond to consideration of only those solutions in 
which the shape of the system does not change (homothetic motions [25]). We thank Sean Gryb for 
drawing our attention to [24] 

28 Laplace presented the principle without a name as more or less obvious. The coining ‘principle of 
indifference’ is due to the economist John Maynard Keynes in his Treatise on Probability (1921). 
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related to the matter content of the universe. However, what we say about prediction 
of mid-point data is in effect closely analogous to one of the most important insights 
in statistical mechanics. Suppose we have a dynamical system that has a known fixed 
energy and consists of a gas of a very large number of ‘elastic balls’ in a box of fixed 
size with walls off which the balls bounce elastically. Then the overwhelming majority 
of microstates, coarse grained and counted using the Liouvillc measure, correspond to a 
very good approximation to a Maxwell-Boltzmann distribution. If one now makes the 
Laplace assumption that all microstates are equally probable, this corresponds to an 
extremely strong prediction of what will be found on a single opening of the box. 

That one does in the overwhelming majority of cases find a Maxwell-Boltzmann 
distribution is explained by the equilibration processes known (and very largely under¬ 
stood) to take place in such a box in which the system is initially in disequilibrium. 
An immense body of experimental data coupled with dynamical modelling therefore 
lies behind the successful prediction just described. Our situation is quite different: 
the single opening of the box corresponds to training our telescopes on the sky to see 
what our universe is like. But this is a one-off test. No cosmologist has made repeated 
observations of universes evolving, let alone equilibrating. 

Thus, we have no accumulated empirical justification for the principle of indifference. 
It is rather a case of necessity being the mother of invention. It forces us to make a clean 
and simple conjecture of the kind that Popper advocated [26]. The sine qua non for us 
to be able to do this is mathematical consistency. This we have established through our 
demonstration that there does exist a well-motivated finite measure for the universe, 
treated as a scale-invariant system of evolving shapes. 

If the nature of our statistical assumption matches standard practice, the use we 
hope to make of it is very different. The fundamental postulate of statistical mechanics 
tells us what to expect at the end of a process: an equilibated state of heat death. Our 
use of the principle tells us what to expect on the Janus surface J and thus provides the 
initial data for the evolution in both directions away from J. It does not tell us how 
things end but how they begin. 


5 Choice of the Primary State Function 


Our system is characterized by a number of state functions (functions only of the current 
state of the system which do not depend on how the system got there). There can be 
up to 6N — 14 such independent functions. In the first instance, we can consider the 
universe as a simple system in the sense of Lieb and Yngvason [27], that is, a system 
that is characterized by a ‘primary’ state function (typically the energy in the examples 
given in [27]). Our primary state function, however, cannot be the energy, because it is 
fixed to zero. Moreover the energy is a dimensionful quantity, and our whole approach 
is based on considering as physically distinct only solutions that differ in shape space. 
This implies that our primary state function has to be dimensionless (scale-invariant). It 
is then natural to take the shape complexity function we identified in [2] as our primary 
state function: 
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(27) 


where f rrris was defined in (2). 

In [2], in fact, we observed that C s is a scale-invariant property of the system that 
characterizes its state rather well: it expresses how clustered/inhomogeneous its instan- 
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taneous shape is. In [2] we proved that C s is bound to grow on average away from 
the Janus point (although it fluctuates, one can bound C s from below and above with 
linear functions of Newtonian time). We also made the observation that C s measures 
well the most striking characteristic of the IV-body solutions, that is, the fact that the 
system becomes more clustered as it evolves away from j. and therefore C s would be the 
simplest way to describe qualitatively the difference between states of a solution near 
to and far from its Janus point. In this paper we make this quantitative and precise. 

In fact, in the first place we use C s to characterize the typicality of mid-point data 
on J. We believe it is the natural, uniquely suitable counterpart of the energy in conven¬ 
tional statistical mechanics. For the inherently nonequilibrium and unconfined relational 
universe, there are no readily available counterparts of the ‘work functions’ (pressure, 
volume and temperature), which played such an important role in creating the frame¬ 
work in which thermodynamics was discovered. It is however noteworthy that the rela¬ 
tional principles that underlie our approach dictate quotienting by scale, which, together 
with dynamical similarity, leads directly to the finite induced measure on PT* S. This 
is the condition that in standard statistical thermodynamics is achieved by restriction 
of the system to a finite volume. The relational universe ‘comes with its own box’. 

This key fact, together with the existence of the primary state function, means that 
an adequate statistical framework does exist in the cosmological context. When we 
discussed Laplace’s principle of indifference, we commented that a known law gives no 
guidance which solution of it will be realized and, for this reason, there appeared to 
be no alternative to ascribing equal probability of occurrence to each of them. It is 
therefore very striking that, given our need to employ Laplace’s principle if we are to 
establish the typicality of solutions, the relational fV-body problem gives us precisely the 
structures we need to do that: the Janus surface J, whose quotient by the gauge orbits 
of Sim is symplectic to T* S; the space of physically distinct solutions, which coincides 
with PT* S; the Liouville measure on T *S and a metric on S which allow us to induce a 
finite measure on PT* S; finally, a uniquely defined primary state function. 

Further state functions can be used to achieve a more restrictive definition of classes 
of microstates considered to be macroscopically indistinguishable. Complexity presents 
itself as the most immediately obvious choice by its simplicity and close analogy with 
energy. 29 There are many ‘finer’ measures of macroscopic distinguishability, e.g., the 
complexity of subsystems or the power spectrum of density fluctuations. In principle, 
one could introduce a complete set of independent state functions, at which point all 
‘entropic ignorance’ is eliminated and one is dealing with individual microstates [7]. 
What one needs are a few state functions that capture the most salient features of 
the evolving microstates. In this paper, we shall merely show what can be done with 
the primary state function: it is adequate to study the typicality of solutions and an 
excellent first diagnostic of the effect of evolution off the Janus surface. 

29 In fact, as we noted in [15, 2], the negative of the complexity is the shape potential , i.e., the part of 
the Newton potential that changes the shape of the universe but not its size. This is exactly what one 
requires of a theory in which overall scale is gauge. Moreover, since the total energy of the relational IV- 
body problem is zero, specification of the complexity, and thus the shape potential, is closely analogous 
to specifying the energy in conventional statistical mechanics. 
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6 Typicality of Solutions and Microstates 

The most important fact that we have so far established is this: every point in PT* S 
uniquely and non-redundantly defines a solution. Using our measure on PT* S, we can 
‘count’ the solutions whose solution-determining mid-point j have given values c of the 
primary state function C s . Consider the natural projection n of points j in J to PT* S 
(this projection associates each point on the surface J with the corresponding gauge 
orbit in T* S, and then projects each cotangent vector in T* S onto the ray in PT* S that 
is parallel to it). Then we define 

£ S oi(c) = Vol {7r (j ) e PT* S / C S G) = c} , (28) 

where the volume Vol is calculated through the measure on PT* S we defined in the 
previous section. As a count of states with given value of a phase-space function, E aol 
is clearly an entropy-type quantity. To avoid confusion, we call it the solution entaxy: 
each j determines a solution (sol); ‘taxis’ is Greek for order and ‘en’ means ‘towards’. 
As we shall see, E sol is a measure of creation by gravity of order out of disorder. 

We use the simple volume rather than a p % log pi integral because we assume each 
possible solution is equally likely to be realized. This is the fundamental postulate 
of statistical mechanics (Laplace’s principle of indifference) applied to solutions rather 
than states. It is also an example of the maximum-entropy principle used to avoid 
bias [7] in a situation in which one has only partial information. In such a case, one 
maximizes the relevant pi log pi expression subject to constraints supplied by known 
prior facts. In our case there is no prior knowledge - given a law, one has no information 
which of its solutions will be realized. It is well known that in the absence of prior 
information the Pi log p.i expression is maximized when all the probablities Pi are equal. 
The pi log pi expression then reduces to the logarithm of the measure. It is important to 
retain the logarithm in the statistical-mechanical definition of entropy because then the 
entropy of independent systems is additive as is found empirically in phenomenological 
thermodynamics. However, in our case we are considering a unique system - the universe 
- so the need for the logarithm falls away. 

Our other important assumption is the measure. Since the whole theory is expressed 
through a phase space with a dynamically distinguished Liouville measure, it seems 
unnatural to choose any other measure. However, the ultimate justification of both our 
assumptions, which are questioned in [19], must be a posteriori through observations, 
which we discuss below. Note, however, that our approach is not subject to the most 
important objection raised in [19] against typicality arguments of the kind we adopt: 
our measure is finite, not infinite, so our probability estimates are unambigous. 30 

We now consider a very important relationship, namely the dependence of the so¬ 
lution entaxy £ sol on the value of the primary state function C s at the Janus point. In 
accordance with our adoption of Laplace’s principle and the induced Liouville measure, 
every shape is equally likely at J. Figure 3 shows clearly that the area on the 3-body 
shape sphere corresponding to shapes near dart-like triangles (with two particles rela¬ 
tively much closer to each other than to the third) is small. In contrast, the moderately 
scalene triangles occupy the bulk of the shape-sphere area. Moreover, most of these 

30 An analogous argument based on the use of a Janus point, or ‘bounce’, is made in [24], However, 
whereas ‘bounce’ is clearly an appropriate word to use in the context of loop quantum cosmology, 
in which a quantum bounce ‘intervenes’ to halt and reverse the collapse inherent in classical general 
relativity, all of the Janus points we consider correspond to smooth unhindered passage through points 
whose existence is mandated by the underlying time-symmetric classical equations of motion. 
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Figure 3: Shape complexity (27) as elevation plot over the equal-mass 3-body shape space. 
The minimum is at the most ‘uniform’ shape (the equilateral triangle) and grows unboundedly 
in ‘cluster formations’ ( r a b <C r ac ,rb c ). The lengths of the equal-Cs curves measure the typi¬ 
cality of the corresponding shapes. For large N, minimal C s correspond to extremely uniform 
(super-Poissonian) distributions [28], large Cs to strongly clustered distributions. 


triangles have C s values not much higher than the minimum at the equilateral triangle. 
This fact is clearly seen in the histogram in Fig. 4. The effect becomes more and more 
pronounced as N gets larger (Fig. 5). 

Thus, on the the basis of the two assumptions we have made, it follows that the 
overwhelming majority of solutions will have values of their primary state function C s 
very near but not exactly at its absolute minimum. Thus, the most typical states of the 
solutions at their respective Janus points will have very uniform spatial distributions of 
the particles. 

The overwhelming majority of solutions will go from C s values very close to (but 
not exactly at) the minimum at J to infinite limits as t —> ±oo. In the continuum- 
approximating limit N —>■ oo, the disorder-to-order effect will be very strong. 

We now explain why £ sol measures creation of order out of disorder. For this we 
must first define the current (instantaneous) entaxy of a solution. So far, we have only 
defined the solution entaxy £ aol . At each point s on a solution, the current entaxy £(s) 
is defined simply as the entaxy corresponding to the shape s: 

£(s) = Vol Mj) e PT *S / CsG) = C's(s)} . (29) 

We noted earlier that solution-determining data away from the Janus point include not 
only the instantaneous shape s and direction d in S but also some quantity that measures 
the ‘bending’ of the solution curve in S at the considered instant. This implies that, for 
each shape s and direction d in S, there exists a one-parameter family of solutions, and 
just one solution of the family has its Janus point at s. Our definition of the current 
entaxy £(s) corresponds therefore to the £ sol of the solution that has its Janus point at 
the shape s. 
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Figure 4: The fraction of the volume of the equal-mass 3-body shape space occupied by config¬ 
urations with different values of the complexity C. The above diagram should be understood 
in the sense of a probability density function: taking an interval C E (Cj, C 2 ), the area below 
the curve Vol (C) between Cj and C -2 gives the fraction of the volume of S occupied by states 
with C between Cj and CV This diagram has been obtained numerically, by sprinkling points 
on S with a probability given by the natural measure induced by Q on S. Notice that the 
curve is, of course, zero for C E (0, C min ). Then it jumps very quickly to a high value followed 
by an initially steep but then much more gradual decline. 


Vol(C) 
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Figure 5: The fraction of the volume of the IV-body shape space occupied by configurations 
with different values of C /C min for five values of N and for a Gaussian distribution of equal- 
mass particles. With this normalization, it is apparent how the distribution gets peaked around 
smaller values of C/C min as N grows, and how the peak becomes more pronounced. 
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We showed in [2] that in every TV-body solution with E = L = 0 the shape complexity 
C s fluctuates around its minimal value in the neighbourhood of the Janus point (where 
D — 0) and then grows with fluctuations but within monotonically rising bounds in 
both directions away from J as t —y ±oo (see Fig. 6). This disorder-to-order effect is 
already present in the 3-body problem, in which more or less chaotic 3-body dynamics 
near J evolves in both time directions into non-chaotic motion of an inertial particle 
and a Kepler pair moving in the opposite direction to the single particle. For TV > 3, 
the system breaks up into subsystems whose centres of mass separate in Hubble-type 
motion as in the inertial system considered in Sec. 2. In the subsystems, energy and 
angular momentum are increasingly well defined and conserved. Overall, the complete 
TV-body system becomes increasingly well ordered. Using the emergent Kepler pairs as 
rods and clocks, internal observers either side of J will find themselves in an expanding 
universe, as illustrated by Fig. 6. For details, see [15, 2], 

Thus in every solution the primary state function C s has minimal values near J 
and then fluctuates between monotonically rising bounds as t —y ±oo (Fig. 6). Since 
low C s values correspond to unstructured uniformly random states and large values to 
well separated subsystems with well-defined conserved quantities, it is clear that, with 
increasing distance from J, gravity inevitably - in all solutions - creates order (and with 
it an arrow of time) out of random disorder. Our new theory using statistics of mid¬ 
point data now suggests we can predict how pronounced the disorder-to-order effect will 
be. 

In Sec. 2, we questioned the possibility of applying standard entropic arguments to 
a universe that is unconfined and therefore can expand freely. This led us to cast doubt 
on the claim that the entropy of the universe is increasing. We have now shown that, 
at least in the TV-body problem, entaxy is a sensible substitute for entropy and that 



Figure 6: Top: plot of the shape complexity C s found by numerical simulation for TV = 1000 
equal-mass particles with Gaussian initial data satisfying E cin = L = P = 0. Bottom: ‘artistic 
impression’, i.e., not found by numerical simulation, of three (scaled) configurations at different 
Newtonian times t corresponding to three points in the solution at the top. 
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it decreases in all solutions away from the Janus point. Any observer must be on one 
or the other side of the Janus point and find herself in a universe that is expanding 
in the direction of decreasing entaxy. In Sec. 7, we will consider how this conclusion 
can be reconciled with the second law of thermodynamics. Before that, we will show 
how one can extract meaningful predictions from our assumptions and some statistical 
considerations. 

Our theory is predictive 

We hold that falsification of our theory of the statistics of mid-point data appears to be 
possible through the predictions it makes for the asymptotic behaviour of our model. 
Consider first the 3-body problem (which is planar if L = 0), whose typical solutions 
we describe in Fig. 7 and caption. Except for zero-measure data at J, there are always 
a Kepler pair and an escaping inertial particle in both t —y ±oo limits. In these limits, 




Figure 7: Left: A typical solution of the 3-body problem in a centre-of-mass inertial frame in 
the case E = L = P = 0 (when the motion is always planar). With the time direction show by 
the arrows, one particle of the pair coming from the right forms a new pair with the particle 
coming from the left. The Janus point is in the middle of the region of non-trivial 3-body 
interaction. Right: the same solution plotted on shape space S. The red and the yellow part 
of the orbit are distinguished by belonging, respectively, to the (nominal) past and the future 
of the point at which the dilatational momentum D is zero and the centre-of-mass moment of 
inertia 7 cm is at its minimum. This point is shown as © at the back of S, where the orbit is 
dashed. Whereas the complexity level-lines are shown for the spiralling yellow half, this cannot 
be done for analogous the red (dotted) half at the ‘back’ of the shape sphere. In this diagram, 
we see clearly the dramatic effect of the dynamical attractors that act on the shape degrees 
of freedom when the ‘veiling’ effect of scale is factored out of the conserved Liouville measure 
on the Newtonian phase space. As with inertial motion, the evolution curves are ‘sucked’ into 
ever-decreasing regions of shape space. However, at the same time, structure is created (in the 
form of a Kepler pair in this simplest example of the 3-body problem). Thus, already in the toy 
universe we also see that Hubble-type expansion is an inevitable attractor in a gravitational 
theory formulated on an extended phase space with scale as a degree of freedom. Moreover, 
the attractor effect is accompanied by ever more perfect formation of structure. A similar 
diagram, for a different solution, can be seen in [15]. 
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two distinguishing scale-invariant quantities stabilize: the eccentricity e of the pair, and 
the angle 9 between the semimajor axis of the orbit of the pair and the coordinate vector 
of the isolated particle relative to the centre of mass of the 3-body system. Numerical 
calculations for mid-point data taken randomly on the Janus surface give characteristic 
distributions of the asymptotic values of e and 9 , as shown in Fig. 8. 



Figure 8: Distribution of the asymptotic values of e and | cos0| for numerical simulations of 
the equal-mass 3-body problem obtained by choosing mid-point data at random on the projec- 
tivized cotangent bundle PT* S with uniform distribution according to the natural measure on 
PT* S; see Appendix A). Our code performs a simulation starting from the Janus point until 
sufficiently late times that a Kepler pair is likely to have formed. Then it identifies the pair 
by proximity (the threshold is that rf 2 < 10 _4 r^ 3 , a = 1,2), and it calculates the eccentricity 

through the formula e = yj 1 + ^^K pair ||L pair || 2 , where E pail and L pair are the energy and the 
angular momentum of the pair, and m red is the reduced mass of the pair. At the same time, 
our code calculates the cosine of the angle between the direction of the isolated particle and 
the semimajor axis of the pair. The latter is determined by looking for local maxima of the 
pair length r\ 2 - 


We can go further and consider the correlations of the asymptotic values of e and 9 
with the value of C s at the Janus point (which is related to the solution entaxy £ aol ). In 
the case of the eccentricity e, we obtain the plot in Fig. 9. 



Cs 


Figure 9: Correlation between the asymptotic values of e and the Janus-point values of C s . 
This plot is the result of ~ 5000 simulations. 
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There is clearly some correlation: if we make a histogram of occurrences of values of 
e and C s ‘cut out’ vertically a some small interval of values of C s , we obtain a different 
distribution if the interval is taken around C s ~ 1.1 C in , n or around C s ~ 1.6C min . This 
is more clearly seen in the case of 9: we can plot the histogram of occurrences of 9 at 
values of C s smaller or larger than 1.5C min , and we get two very different distributions, 
as shown in Fig. 10. 
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Figure 10: Distribution of the asymptotic values of |cos0| for solutions with C s smaller or 
larger than l.hC^in at the Janus surface J. 

At small Janus-point-complexities, we get a distribution with a peak at 9 = 0 and 
a plateau for any other value of 9. For large complexities, on the other hand, the small 
values of 9 are penalized in favour of 9 — 0. This corresponds to an, as yet admittedly 
weak, prediction: since small values of C s at the Janus point are statistically favoured, 
our model predicts a distribution of 9 more like the one on the left in Fig. 10. If, in 
our hypothetical 3-body universe, values of 9 far from 0 were measured, that would 
correspond to evidence against our model. 

Of course the predictions that one can make with the 3-body problem are extremely 
limited and weak, and serve mainly as a proof of principle. 

However we conjecture that this dependence of the distribution of asymptotic data 
on E sol becomes more marked for some shape variables in the IV-body problem for large 
N, for which there will be vastly more scale-invariant quantities that can be correlated 
with the solution entaxy. If it is correct, then the hypothesis that the solutions with high 
E sol are the ones most likely to be realized can be tested on the basis of the asymptotic 
behaviour. 31 


7 Entropy and Emergence of the Second Law 

Our discussion of typicality of the shape-dynamical description of a gravitationally dom¬ 
inated universe raises some questions in connection with the second law of thermody- 

31 In the IV-body case, the system will in general break up into several clusters and one could, for 
example, make statistical studies of how the average asymptotic quantities (say 1+3 as opposed to 2+2 
decomposition for N = 4) depend on the typicality of the data at J. For large IV, predictions could 
be made for many observable quantities for typical and atypical mid-point data. Observers in one or 
the other of the asymptotic limits (as we have noted, they cannot be in both) will always observe an 
arrow of time and, if our conjecture is correct, be able to establish the extent to which the universe in 
which they find themselves is typical or atypical. Any IV-body universe that is maximally typical at J 
will, by virtue of being so, evolve away from J to states which are highly atypical and carry information 
confirming the maximal typicality at J. This is exactly opposite to the behaviour mandated by the 
second law of thermodynamics, according to which atypical states evolve to typical ones. Obviously, 
we must explain the discrepancy, to which we now turn. 
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namics. First, entropy is generally understood as a measure of typicality and the second 
law of thermodynamics states that entropy never decreases. In Sec. 2 we already made 
what we think is a convincing case for the inadequacy of thermodynamical entropy as 
a measure of typicality of the state of the whole universe. Typicality on shape space is 
different from thermodynamical entropy, which depends on dimensionful quantities and 
can only belong to subsystems of the universe. This is well illustrated by the exam¬ 
ple of free expansion of an ideal gas (inertial motion) shown in Sec. 2: evolution does 
not correspond to an increase in typicality, because as the gas expands correlations are 
established between positions and momenta that are not present at the Janus point. 32 

In the rest of the present paper up to this point, we showed how our entaxy captures a 
suitable notion of typicality of the state of the universe and how its tendency to decrease 
is linked to the gravitational arrow of time we identified in [2]: it is a consequence of 
the presence of dynamical attractors in shape space (see Fig. 7). We believe that this 
resolves the tension between typicality and the arrow of time [13, 19, 6], and also opens 
the way for ‘retrodiction’ of the likelihood of past conditions given a knowledge of the 
present ones. This was deemed impossible, for example, by Shiffrin and Wald [19], in 
a universe governed by the second law. The example given in Sec. 6 shows this is in 
principle possible. 

What we have obtained so far is in not in tension with the standard formulations of 
the second law of thermodynamics. On the contrary, the growth of thermodynamical 
entropy will be seen to emerge as a characteristic pattern in all solutions of our model. In 
fact, we now want to examine the behaviour of the subsystems as they form progressively 
and become ever better isolated from the rest of the universe. We will argue that one can 
assign to them an entropy, defined in the conventional way, which is initially low and then 
increases. Thus, we have two different kinds of entropy: one for the universe (entaxy) 
and one for the subsystems that form within it: while the former decreases, the latter 
increases, matching the standard formulation of the second law of thermodynamics. 

Entropy in self-gravitating systems 

The notion of entropy we refer to is the microcanonical entropy , which measures our ig¬ 
norance of the microscopic state of a system. It is mathematically realized by the Boltz¬ 
mann relation S = k log W, where W, in the formulation due to Planck, is the volume 
of phase space occupied by the states that are compatible with some prescribed values 
of the state functions (the energy above all). This is, according to Gross [29], still today 
“the deepest, most fundamental, and simplest definition of entropy”. Microcanonical 
entropy is suitable for the description of systems with phase transitions or non-extensive 
systems with long-range interactions, for which the canonical Boltzmann-Gibbs statis¬ 
tics fails 33 [30] and the thermodynamical limit is often singular [29] . In particular, the 
microcanonical entropy is the only option to describe self-gravitating systems [31, 30], 

32 Moreover, if adiabatically expanding walls which allow the gas to continually re-equilibrate are not 
present, there is no useful definition of entropy (already Gibbs’ entropy is ill-defined). 

33 The specific heat of any system described by a canonical ensemble will be positive-definite, c v = 
f3 2 ((E — ( E )) 2 ), therefore this ensemble cannot describe systems with negative specific heat. More 
physically, the canonical ensemble is supposed to describe a system in contact with a thermal bath. If 
the system has negative specific heat, it will not come into equilibrium with the bath because in such 
a situation heat flows from hot bodies to cold ones, in violation of Clausius’ formulation of the second 
law. 
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because they have a negative specific heat. 34 

Microcanonical entropy is not well defined in Newtonian self-gravitating systems 
unless both a UV and an IR cutoff are imposed to ensure that the hypersurface of 
constant energy has a finite volume. In fact one needs an IR cutoff, for example a 
reflecting wall surrounding the system; otherwise the configuration part of the constant- 
energy surface is non-compact (one can always put a particle arbitrarily far away with 
a small enough kinetic energy in such a way that the total energy does not change). 
Also necessary is a UV cutoff preventing Newton’s potential from diverging when two 
particles come arbitrarily close to each other. This can be realized either with a repulsive 
hard core or with a flattening of the potential below a certain distance. Either way, the 
purpose of this cutoff is to introduce a minimum for the potential energy, which, in turn, 
implies a maximum for the kinetic energy on the constant-energy hypersurface. This 
makes the hypersurface compact. 

Now, it is clear that the nature of our toy model of the universe would be deeply 
changed by putting it inside a box with reflecting walls. Therefore an IR cutoff cannot be 
imposed and microcanonical entropy cannot be defined for the whole universe. On the 
other hand our model spontaneously produces bound subsystems which are effectively 
‘boxed’ by their own gravitational potential well for a certain time until a particle 
escapes, leaving a smaller system that, in its turn, is effectively boxed. For these systems, 
once they differentiate from the rest of the universe, it makes sense to speak about a 
microcanonical entropy and an IR cutoff. Then, for these subsystems, all the well- 
established results on the thermodynamics of self-gravitating systems apply [30], and 
entropy growth is observed as a characteristic tendency of the subsystems to evolve 
towards local metastable maxima of the entropy that can escape the UV problem for a 
considerable time. This fragmentation into subsystems for which the second law holds 
is an example of the formation of ‘branch systems’, proposed by Reichenbach [33, 6] 
to solve the puzzle of time-reversal. 35 Davies ([6], p. 69) describes branch systems as 
follows: 

Branch systems are regions of the world which separate off from the main envi¬ 
ronment and exist thereafter as quasi-isolated systems, and usually merge once 
again with the wider environment after a sufficient time. Examples of this sort 
are countless, but one will suffice to remove any misunderstanding. When we take 
an ice cube, and add it to a lukewarm drink, and watch the ice cube melt, the 
system ice+drink only comes into being after the event. It simply did not exist 
as a quasi-closed system beforehand. Also it will be seen that it really is a quasi- 
closed system in the sense that the melting process can be perfectly adequately 
described without recourse to interaction with the outside universe. 

For avoidance of any possible misunderstanding, we note that the human intervention in 
Davies’ example is immaterial; an iceberg falling off Greenland into the Atlantic would 

34 A negative specific heat is a common feature of all long-range potentials V ~ r~ a , a < 2. This can 
be proven assuming that the system under consideration is (meta-)stably bound, so that the dilatational 
momentum is constant on average [32]. Then the virial theorem holds, 2(T k ) ~ — cx(V), and one can 
write (E) ~ (1 — ^)(T k ), which gives a negative specific heat c v = Jyfjy if a < 2. 

35 Despite the utility of Reichenbach’s notion, especially as characterized by Davies, we think Re¬ 
ichenbach completely fails to explain the origin of the arrow of time. He assumes that the universe 
is subject to Poincare recurrences and at no stage even mentions gravity let alone attributes to it a 
significant role. Davies does not commit himself to the success of Reichenbach’s attempt, merely saying 
“Reichenbach asserts”. 
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serve just as well. 


36 

A dynamical derivation of the second law 

We can now describe our task. We have to explain three facts. First, throughout the 
universe we see evidence for copious formation of branch systems as described by Davies. 
Second, the overwhelming majority of them evolve in essentially the same characteristic 
way that can be called increase of entropy; in each such branch system, this growth 
defines an arrow of time. Third, all the arrows associated with the formation of such 
branch systems point in the same direction. These are the facts we must explain, at 
least in a reasonably realistic toy model. 

We will of course attempt to do this by means of systems with Janus-point solutions. 
Now there are numerous systems described by equations that are time-reversal symmet¬ 
ric and have Janus points, but three we mentioned earlier - inertial motion, water waves 
and systems with repulsive gravity - do not form branch systems. In contrast, attractive 
gravity does. 37 

Our immediate aim is to show that a typical universe is perfectly compatible with 
the formation of branch systems with arrows of time that, on each side of the universe’s 
Janus point, all point in the same temporal direction. This we can do using the classical 
iV-body problem. 

We will summarize here in some more detail this process of fragmentation and how 
it gives rise to subsystems for which the second law holds. 

• According to corollary 1 of [34], a typical solution iV-body solution at late times 
(from our perspective, away from the Janus point in either time direction) will 
fragment into ‘subsystems’ whose centre of mass drifts away from R cm linearly in 
Newtonian time t. 38 Each subsystem consists of several ‘clusters’ which drift away 

36 Albert [3] is critical of the branch-system concept (p. 88/89), seeing two main problems in it. First, 
how is it possible “to decide at exactly what moment it was that the glass of water with ice in it first 
came into being?” Second, he sees arbitrariness not only in the temporal origin of a branch system 
but also in its spatial extent. Why not include the room that contains the glass of water with the ice, 
the building containing the room and so ad infinitum ? In fact, a few more words by Davies could have 
answered these criticisms. First, the moment the ice enters the water is, macroscopically, a well-defined 
instant that creates locally a manifestly strong departure from equilibrium. Second, one can, with 
Albert, identify (define) infinitely many branch systems larger in both the temporal and spatial extents 
that all include the glass, water and ice. These will have boundaries, which need not be physical as 
long as they are well defined, in both time and space. Now because the number of surface degrees of 
freedom scale as the square of the linear dimensions while the number within a volume scale as the cube 
of those dimensions, the positions of the temporal and spatial boundaries have no sensible effect on the 
melting of the ice. Whatever else may be happening in these extended branch systems, the melting of 
the ice in Davies’ minimal one will contribute an entropy increase in all of the extended ones. 

37 Attractive gravity is crucial and its role intriguing. One can barely imagine that the universe is 
like steam in a cylinder ■ it is surely unboxed whatever law governs its behaviour. If we are to obtain 
a finite measure and hence formulate meaningful statistical statements of the kind made possible by 
the containers always presupposed in historical studies of thermodynamics and statistical mechanics, 
we cannot avoid quotienting by scale. This move is supported by the fact that within the universe 
all observations are ultimately of ratios. However, once branch systems, above all Kepler pairs, form 
through the action of gravity, dimensional units are defined by them to an ever better accuracy. More¬ 
over, and this completes the circle, the gravitational potential of each negative-energy branch system 
ensures that it remains effectively confined to a finite region of space - it creates its own box. Note 
that if and when a particle does escape from the box, it carries away positive energy. This leads to the 
formation of a new branch system that is even more tightly bound. 

38 See [1, 15] for a definition of Newtonian time purely in intrinsic shape-dynamical terms. 
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from each other at most as f 2 / 3 [34], 

• The fact that the system fragments into clusters has the important consequence 
that the rest of the universe provides units to describe the state of a cluster. To 
see this, we can, for example, take an effectively isolated Kepler pair as a reference 
system. We use its mass to define the reference solar mass and use its semimajor 
axis as the reference astronomical unit and its orbital period as the reference year. 
Moreover, we can use its center of mass to define a coordinate origin and define 
a frame by the orbital plane and the semimajor axis of the Kepler pair. It is 
thus possible to discuss clusters of the iV-body system using units of length, time 
and energy. See [1] for a more detailed description of the procedure for defining 
reference frames, a size and a time standard from relational/shape data. It is 
these definitions that enable us to describe clusters in terms of size, separation, 
linear momentum, angular momentum and energy. With increasing distance from 
the Janus point, each cluster will therefore be characterized by physical quantities 
that become progressively better defined and conserved [15]: energy, linear and 
angular momentum. 39 These conserved quantities can be used as state functions 
for a statistical-mechanical description in terms of the microcanonical ensemble 
[29], 

• The subsystems virialize, as their moment of inertia grows at most like t 4 / 3 , 
and a necessary and sufficient condition for the strong form of the virial theorem 
described in [32] is that limi -2 1 = 0, which is clearly satisfied. 40 Then the second 
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law will be seen to hold, in the sense that coarse-grained descriptions of the system 
(given, for example, by the mean-held approximation 41 ) will assign an increasing 
Boltzmann entropy to the vast majority of physical solutions of the system [30]. 

• The statistical mechanics of self-gravitating systems is extremely rich and is closely 
related to that of systems with first-order phase transitions. The system admits 
two distinct phases: a hot diffuse gas which can effectively be considered an ideal 
gas with no potential energy, and a collapsed core with very little kinetic energy. 
These two phases coexist in actual solutions and the process of evaporating parti¬ 
cles from the core into the diffuse component, which has the effect of collapsing it 
further, is favoured because it is in the direction of growing entropy [30]. This is 
the basis of the so-called gravo-thermal catastrophe [36]. The configurations that 

39 These asymptotically conserved quantities are a consequence of the existence of asymptotic Eu¬ 
clidean symmetries for clusters (when a cluster is infinitely far from the others it acquires an effective 
Lagrangian that is invariant under rotations and space and time translations). Noether’s first theorem 
then applies with ever better accuracy. 

40 In contrast, the universe has lim t~ 2 I = const > 0, which is why we cannot treat it the same way 
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as its clusters. 

41 The mean-field approximation only considers the one-particle distribution function /(p, q) and 
ignores all correlations between particles (n-point functions). It can be shown to arise as the dominant 
contribution to the microcanonical entropy from a saddle point approximation. The entropy S B = 
f d 3 pd 3 g / log / associated with this coarse graining is the Boltzmann entropy; it shows a clear tendency 
to grow in the vast majority of solutions, which, starting from generic initial conditions, will evolve 
towards configurations which maximize S B [30]. This is the sense in which virialized self-gravitating 
systems satisfy the second law. The complications due to the long-range nature of gravity relate 
to the non-existence of entropy-maximizing configurations: there are only local maxima (metastable 
equilibria) and ultimately no global maximum of S B . As Sanmartin [35] puts it “a runaway in entropy 
is a common fate of gravitationally bound systems”. 
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realize local maxima of the entropy for a system with given energy and angular 
momentum possess a rich and interesting phase diagram [31]. 

This shows how the overwhelming majority of branch systems formed in our toy 
model evolve in essentially the same characteristic way that can be called increase of 
entropy. 

ft remains to explain why the various thermodynamic arrows of time present in 
subsystems of the universe, the existence of which constitutes the experimental basis of 
the second law, all point in the same direction as the global gravitational arrow. How 
does it come about that, in our current epoch, we observe such a plethora of low-entropy 
branch systems in all of which the entropy is increasing in the same direction? How do 
all these unidirectional low-entropy branch systems arise? 

The only widely accepted answer to this question is the past hypothesis 42 - that the 
universe began in the big bang in a state of extraordinarily low entropy. Our theory of 
mid-point data in the self-gravitating iV-body problem suggests that the problem is an 
artefact of the transfer of a conceptual framework suitable for the study of steam engines, 
in which the steam is confined and its self-gravity plays no role, to the self-gravitating 
universe. 

In our framework, the perceived problem with time-reversal symmetry that leads to 
the ‘second-law paradoxes’ is eliminated. Because observation is necessarily restricted to 
part of the universe that is either on one or the other side of its Janus point, the effective 
law which governs the observed half is not time-reversal symmetric but asymmetric. 
Coupled with the fact that gravity is attractive and tends to clump matter, the formation 
of branch systems that become ever better decoupled from the rest of the universe is 
a dynamical necessity. The only role that statistics plays in this process relates to the 
specific parameters of the branch systems. The fact that they form and with time are 
described ever better by Hamiltonian equations, so that (negative) energy and centre- 
of-mass angular momentum are conserved more and more accurately, is inescapable. 

It now only remains to show that the all-pervasive dynamical arrow of time thus 
imposed will ensure that the overwhelming majority of the created branch systems 
will begin life with a relatively low Boltzmann entropy that will then increase with 
increasing distance from the Janus point. Since realistic self-gravitating systems cannot 
have eternally stable maximum-entropy states, we merely note that in our scenario 
virtually all branch systems consisting of enough particles to make statistical arguments 
justified will form with some entropy and then evolve to states of higher entropy in 
which they exist for a relatively long period. Globular clusters provide a good example 
of what we have in mind: they typically consist of ~ 10 6 stars and are mostly well fitted 
by the Michie-King model [37] with a core-halo structure. Many of these were formed 
very early in the universe and have been slowly evaporating ever since. 

There is also the issue of why an initially low entropy virtually never decreases as 
it would if all the velocities of an entropy-increasing branch system were reversed at 
the time of its formation. Statistically, both velocity directions are equally possible. 
However, it is now well understood in conventional statistical mechanics [6] that, in the 
overwhelming bulk of cases, an initial decrease of entropy will almost immediately be 
transformed into an increase and, moreover, so rapidly as to be effectively unobservable. 
The reason for this is that the number of accessible microstates decreases very rapidly 
with the entropy, so that an entropy-decreasing system very soon ‘finds its way’ for 
purely statistical reasons into a phase-space region of higher entropy. 

42 The coining is due to Albert [3]. 
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As regards the other arrows of time, we have already suggested (in footnote 4), that 
retarded potentials will be observed on either side of a Janus point. Further, it is widely 
accepted that the psychological arrow of time is determined by the thermodynamic 
arrow, so that should be explained if our overall picture is correct. Another important 
arrow of time is the one associated with ‘collapse’ of the quantum wave function. On 
this point, we will merely note that the manner in which the information associated 
with entanglement ‘leaks into the environment’ as a result of decoherence is strikingly 
similar to thermodynamic ‘loss of information’, so that too may have its origin in a 
Janus point. 

To summarize, the big picture is that, on either side of its Janus point, the decreasing 
entaxy of the universe will define an unambiguous arrow of time with which the emergent 
arrows of the overwhelming majority of the branch systems that form will naturally be 
aligned. We have, of course, said nothing about non-gravitational forces, which supply 
much ‘cleaner’ examples of entropy growth than the astrophysical examples we have 
mentioned. However, we would argue that the formation of any branch system that 
then increases its and its environment’s entropy always has its ultimate cause in the 
clustering effect of extaxy-decreasing gravity. It clearly enabled the sun and earth to 
form, which in turn led to the creation of coal, the invention of steam engines, and the 
discovery of the second law. 

8 Conclusions 

In the body of the text, we have said virtually nothing about a key question: is our proof- 
of-principle Y-body model a good indicator of what can be expected in the framework of 
general relativity (GR) or some similar relativistic theory? As we have already written 
about the architectonic similarities between GR and the Y-body problem in [38, 15, 2], 
we will here only make some salient points. 

• The initial-value problem of the relational Y-body problem, in which a point and 
a tangent vector in shape space determine a solution, is closely analogous to that 
in vacuum GR, in which a point and direction in conformal superspace determine 
a solution [39, 38]. It follows that it will be possible to specify scale-invariant 
(conformally invariant) mid-point data if suitable Janus points are present in the 
physically relevant solutions of GR. Then our theory of the typicality of universes 
should be applicable. 

• In fact, Janus points are present in some classes of GR solutions, for example in 
big-bang-big-crunch solutions. However, in those cases they are at maximum ex¬ 
pansion, whereas in the Y-body problem they are at minimum expansion. There¬ 
fore, more promising is what happens at the big-bang. Since overall scale is a 
gauge degree of freedom in a truly relational theory, we do not need to worry 
about the scale singularities that are said to indicate breakdown of classical GR. 
We merely require regular continuation of the conformal 3-geometry through the 
big-bang to a universe ‘on the other side’, in which its time will have the opposite 
direction to ours. This is the most immediately relevant issue to resolve. 

• In GR, the three-dimensional conformally invariant Yamabe constant is a quantity 
closely analogous to our complexity and could therefore play the role of a primary 
state function [15]. 
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• In its broad features, above all in gravitational clustering and the associated for¬ 
mation of branch systems, the history of our universe since it became transparent 
at the surface of last scattering does seem to be well reflected in the asymptotic 
behaviour of the IV-body problem. This is circumstantial evidence in support of 
our proposal. 43 

Besides these points, we see hope that we are on the right track in the sheer natural¬ 
ness of an unconfined system like the universe (governed by a time-reversal symmetric 
law) having a Janus point. Nearly a decade before we realized the significance of the 
one long known to exist in the IV-body problem (since 1772 in the 3-body problem!), 
Chen and Carroll [12] conjectured the existence of such a point. Moreover, after we 
realized that solution-determining data could be specified at such a Janus point, we 
found that this idea had been partially anticipated by Carroll with his ‘middle hypothe¬ 
sis’ [13]: if the underlying law of the universe is time-reversal symmetric, one could well 
expect its solutions to be qualitatively symmetric about a ‘middle’. We also discovered 
that Ashtekar and Sloan (footnote 27) had already noted how solution-determining data 
could be specified at a quantum ‘Janus bounce’, though apparently without recognition 
that this would give a one-past-two-futures scenario. 

In view of these facts, we find it especially striking that already Boltzmann (or 
rather his assistant Dr Schuetz) [5] had, 120 years ago, proposed a one-past-two-futures 
resolution of the arrow-of-time paradox in the context of a universe subject to Poincare 
recurrences and hence necessarily confined within a ‘box’ (since Poincare’s theorem 
only holds if the accessible phase space has a bounded measure). Is it not strange, 
given the long-known ‘Boltzmann-brain’ problems [40] associated with this idea, that 
the implications of a Universe Unbound were never, it seems, contemplated? 
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Appendix A. The measure in the 3-body problem 

In this appendix we will explicitly perform the phase space reduction in the case of 
three bodies. The extended phase space is R 6 , with coordinates iq,r 2 ,r 3 £ R 3 for 
the particle positions, and p 4 ,p 2 ,p 3 € R 3 for the momenta. It is well known that 
if the angular momentum L = y )^ =1 r a x p“ is orthogonal to the plane identified by 
the three particles, L x ((rq — ra) x (r 2 — ra)) = 0, then the three particles never leave 

43 Taken together, these points suggest that key properties of the microwave background, such as the 
temperature equality to ss 10 _5 in regions that could never have been in causal contact, could have an 
acausal explanation in terms of Laplace’s principle applied to the mid-point data of a dynamical system 
with a very large number of degrees of freedom. If the universe does have a Janus point, Fig. 5 suggests 
that its mid-point data will correspond to a spatially very uniform state with however microscopic 
fluctations that can be expected to grow in the initial evolution away from the Janus point.We may 
also mention that all Janus-point solutions of our self-gravitating model will have an at least weakly 
expressed emergent second law. The potential value for cosmology of our theory of typicality is in the 
predictions it will allow for the ‘initial conditions’, i.e., mid-point data, of the universe. 
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that plane during the evolution. In our case the angular momentum L is zero and 
therefore the problem is planar. We can then assume that the position and momenta 
are two-component vectors rq, r 2 , r 3 , p 1 , p 2 , p 3 € R 2 . The degrees of freedom are six, 
but three are gauge, corresponding to the translations and rotations on the plane of 
the motion. We have to restrict to the hypersurface P = Lj_ = 0, where Lj_ is the 
remaining nonzero component of the angular momentum, and then we have to quotient 
by the transformations generated by these constraints. It turns out that in this case we 
can take the ‘royal road’ of explicitly identifying a sufficient number of gauge-invariant 
degrees of freedom (observables), and perform a coordinate transformation in phase 
space that separates them from the gauge degrees of freedom, making them orthogonal 
coordinates. 

To deal with translations, we define the mass-weighted Jacobi coordinates: 


Pi = 


- / mim_ ( r _ r ) 

7711+777,2 V ^ ’ 


0 __ +3 (uii+m;) ( _ mi ri+m 2 r 2 

Y mi+m2+m3 y ^ mi _|_ m2 J > 

P3 = Vmi + m2+m 3 {mi ri + m 2 r 2 + m 3 r 3 ). 
The transformation to them is linear and invertible, 

p a = M a b r b , det M = v / miro 2 m 3 , 
and it is therefore easy to ford a canonical extension for it: 

n a — (M~ 1 ) a b p b — ((M~ 1 ) T ) b a p b . 

The inverse transformation is 

p a = M a b K b = {M T ) b a K b . 

Note that the inverse matrix is 


(M 


— 1 \a 


( J m.\(mi+rn2+rnz) / m\(mi+m2+mz) 
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and has a constant column. It is the column of tc 3 , 
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+ = 


+1 + m 2 + m 3 
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which is therefore proportional to the momentum constraint and vanishes. The coor¬ 
dinates p 3 are gauge; they are the coordinates of the center of mass. The other two 
momenta are 


1 


+ = 


mip — m 2 p 
y/mim 2 (mi + m 2 ) ’ 


(mi + m 2 ) p 3 - m 3 ( p 1 + p 2 ) 
a/( mi + m 2 ) m 3 (r+ + m 2 + m 3 ) 
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The canonical transformation leaves the Poisson brackets invariant: 


(37) 
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The kinetic term is diagonal in the momenta K a , 
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a =1 
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a =1 6,c=l 
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as is the moment of inertia, 
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The inertia tensor also takes a particularly simple form: 


(38) 


(39) 
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a=l 
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We are left with four coordinates pi, P 2 and momenta K l , k 2 , and a single angular 
momentum constraint (the component perpendicular to the plane of the triangle): 




||(n -r 3 ) x (r 2 — r 3 )|| 


(41) 


a =1 


equilateral triangle 



Figure 11: The shape sphere of the equal-mass 3-body problem. Every point on the 
sphere (defined as a constant-||w|| surface) is a triangle. Points at the same longitude 
with opposite latitudes correspond to mirror-conjugated triangles. At the poles (the 
intersections with the W 3 axis) we have the equilateral triangles, while on the equator 
(the red circle, W 3 = 0) we have the collinear configurations. Among them, there are 
six special ones: three binary collisions (red dots, one of which is on the W\ axis), 
and Euler configurations (white dots), in which the gravitational force acting on each 
particle points towards the center of mass and has a magnitude such that, if the system 
is prepared in rest at one of these configurations, it will fall homothetically (without 
changing its shape) to a total collision at the centre of mass. The same thing happens 
at the equilateral triangle (for all values of the masses, as Lagrange showed). 
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where with the vector product between two 2 -dimensional vector we understand a scalar 
a x b = a x b y — a y b x . The coordinates 

wi = \ (IIPill 2 - IIP 2 || J ) , w 2 = pi-p 2 , u >3 = Pi x P 2 (42) 

are invariant under the remaining rotational symmetry and therefore give a complete 
coordinate system on the reduced configuration space. Notice that W 3 changes sign 
under a planar reflection (changing the sign of one of the coordinates, say x, of both 
Pi and p 2 ) while W\ and w 2 remain invariant, and therefore the map W 3 —> — W 3 relates 
triangles conjugate under mirror transformations. This also has the consequence that 
the W 3 = 0 plane contains only collinear configurations (whose mirror image is identical 
to the original, modulo a planar rotation). This has nothing to do with 3D reflections 
(obtained by changing the sign of all components of every Euclidean vector). In fact 
triangles are invariant under such parity transformations, because their parity conjugate 
is related to the original by a non-planar rotation. 

The Euclidean norm of the 3D vector w = (wi,w 2 ,Ws) is proportional to (half) the 
square of the moment of inertia 

INI 2 = \ (llpill 2 + IIP 2 II 2 ) 2 = (43) 

so the angular coordinates in the three-space (wi,w 2 ,ws) coordinatize shape space, 
which has the topology of a sphere [23]. We call it the shape sphere , and in Fig. 11 we 
describe its salient features. 

The norms of the original Jacobi coordinate vectors can be written as 

IIpiII 2 = \Jw\ + w\ + w\ + Wi , I|p 2 || 2 = \Jwl + wl + wl~wi, (44) 

and therefore the full vectors are specified by 


Pi = 
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(45) 


where 5 = arctan ^ is the angle between pi and p 2l and 6 is an overall orientation 
angle which is not rotation-invariant and therefore is not fixed by the specification of 
the coordinates W\,W 2 ,W 3 . Now we want to find the momenta conjugate to w. To do 
so, we consider the symplectic potential 

3 3 

© = p“ • dr a = ^ ( 46 ) 

a= 1 cl— 1 


If we replace p a with their expressions in terms of w from Eq.(45), we get 
0 = z 1 dtci + z 2 dw 2 + z 3 dw 3 + ft 3 ■ dp 3 + L± dd, 
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so we now have a complete canonical transformation from the coordinates (r 1; r 2 , r 3 ; 
p 1 , p 3 , p 3 ) to (u>i, W 2 , W 3 , p 3 , 9; z 1 , z 2 , z 3 , k 3 , L±). The Poisson brackets in these 
coordinates are, as they should be, canonical: 


{z a ,w b } = S a b , {L ± ,e} = 1, {z a ,L ± } = 0, {z a ,z b } = 0, 

{z a , k, 3 } = 0, {9, K 3 } = 0, {z a , /0 3 } = 0, { 9 , = 0. 

In the new coordinates, the kinetic energy decomposes as 
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where the three collision longitudes are [1] 
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If we call 0 the azimuthal and 0 the polar angle on the shape sphere, then = 
cos 0 cos 0 and p|- = sin0 cos 0, and on the constraint surface L± = k 3 = 0 the 
Hamiltonian constraint takes the form 


H = MU \\z\ 


Cs(0,0) 

V¥4 


C s (0,0) = ^ 


a<b 


(m a m b ) 2 (m a + m b ) 2 
a/1 - COS 0 COS(0 - (f>ab) ’ 


(53) 


where C s is the 3-body complexity function. Finally a short calculation reveals that the 
dilatational constraint takes basically the same form in the new coordinates: 


D — W r Q - p" — 2 w ■ z + K? ■ p z . 
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We are now ready to specify our measure on shape phase space. We will gauge-fix 
translations and rotations simply by setting 9 = 0 and p 3 = 0. Moreover the dilatation 
constraint D and the Hamiltonian constraint H will gauge-fix each other: 


a= / \F-P\5 3 (p 3 )6 3 (K 3 )6(9)6(L ± )5(H)5(D)(der 3 , 


where the Faddeev-Popov determinant reduces to 
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while the delta-function of the Hamiltonian constraint can be written as 
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So the product of the two gives 


F-P 6(H) = v/pp (\/M - C‘ /3 ||iir 2/3 ) , (58) 


and we then have 
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If we now perform a last canonical transformation to polar coordinates in shape phase 
space, 

w = w(sin 0 cos 0, sin 0 sin 0, cos 0), 


then 
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which finally leads to 


Notice that we would have got the same result had we used any other gauge-fixing of D 
than H , for example I cm = const. 

Our cr is the natural measure on T* S, the cotangent bundle to the shape sphere. It 
is not integrable, because z^ and z^ do not have a compact domain of definition (or, in 
other words, the cotangent space to a point in shape space is not compact). The last step 
is to recall that solutions with the same mid-point data except for a rigid rescaling of all 
shape momenta (z^,z^) —» const ,(z^,z^) are all similar and therefore indistinguishable 
on shape space. Then there is an equivalence class on T *S of physically indistinguishable 
mid-point data, and we would like to define our measure on the quotient of T *S with 
respect to that equivalence relation, which is PT* S, the projectivized cotangent bundle 
to S. 

Now there is no unique way of inducing a measure on PT* S, but there is a distin¬ 
guished one. Consider the scale-invariant mass metric on the extended configuration 
space: 
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This metric induces a quotient metric on S [1], 


ds 2 = d0 2 + sin 2 0 d0 2 , 
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which, in turn, induces a norm on the associated cotangent bundle T* S: 

||~|| 2 = (,^) 2 + (^) 2 csc 2 0. (66) 

We may then use the space of unit-norm covector fields ( z ^) 2 + {z^) 2 esc 2 if = 1, which 
is the double cover of PT* S, as a redundant description of our solution space. Each 
solution maps to two conjugate points on this space. The conjugacy relation is simply 

So we may dehne our final measure on the space of solutions as 


J 8 ((z^) 2 + (z^) 2 esc 2 "0 — 1) di/jdcfidz^dz^ 

/ 8 ( 'Z 2 — l) sin 0 d^)d(f)dZdx = sin0d0d0dx, 


(67) 


where Z = \J ( z ^) 2 + ( z < t’) 2 esc 2 0 and x — arctan j • This measure is now finite 

since x £ [0, 27r). 
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